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ABSTRACT 

Type II bursts are thought to arise from instabilities in the accretion flow onto a neutron 
star in an X-ray binary. Despite having been known for almost 40 years, no model can yet 
satisfactorily account for all their properties. To shed light on the nature of this phenomenon 
and provide a reference for future theoretical work, we study the entire sample of Rossi X-ray 
Timing Explorer data of type II bursts from the Rapid Burster (MXB 1730-335). We And that 
type n bursts are Eddington-limited in flux, that a larger amount of energy goes in the bursts 
than in the persistent emission, that type II bursts can be as short as 0.130 s, and that the 
distribution of recurrence times drops abruptly below 15-18 s. We highlight the complicated 
feedback between type II bursts and the NS surface thermonuclear explosions known as type I 
bursts, and between type II bursts and the persistent emission. We review a number of models 
for type II bursts. While no model can reproduce all the observed burst properties and explain 
the source uniqueness, models involving a gating role for the magnetic held come closest to 
matching the properties of our sample. The uniqueness of the source may be explained by a 
special combination of magnetic held strength, stellar spin period and alignment between the 
magnetic held and the spin axis. 

Key words: stars; neutron - X-rays: binaries - X-rays: bursts - X-rays: individual: MXB 
1730-335 


1 INTRODUCTION 


X-ray bursts observed in accreting neutron stars (NS) in low-mass 
X-ray binaries (LMXBs) come in two types, sometimes similar in 
appearance (~ 1 min time scales, 10^® erg s“^ peak luminosities) 
but thought to be due to very different mechanisms, because of the 
large difference (by factors ~ 10^) in average power. 

Type I bursts arise from the heating and cooling of the NS 
photosphere after a thermonuclear shell flash of accreted material 
(see reviews by Lewin, van Paradijs & Taam|1993|rStrohmayer &| 


|Bildsten|2006|[^lloway et al.|2008| . This is consistent with ther¬ 

mal spectra of varying temperature and roughly constant emitting 
area. Despite some remaining open questions (e.g., Cornelisse et al.| 


|2003[ jSuleimanov, Poutanen & Wemer|20lT||Watts|2012 l a gen 

eral understanding of the physics underlying type I bursts has been 
reached, also thanks to the wide range of circumstances in which 
they occur (type of LMXB, accretion rates, etc) provided in partic¬ 
ular by the relatively large amount of bursters (over a hundred). 
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Type II bursts remain largely a puzzle. Soon after their discov¬ 
ery, jH^San^^arehajT&Tewj^ljTOT^ showed that their short re¬ 
currence times rule out the possibility that they are powered by nu¬ 
clear burning, ascribing the source instead to gravitational potential 
energy. Little more is understood about their origin, also because 
they are only confidently detected in two sources: the Rapid Burster 
(MXB 1730-335, hereafter RB), a recurring transient which is the 
object of the present study and which also exhibits type I bursts, 
and the Bursting Pulsar (BP), a slowly rotating (2 Hz) accretion- 
powered pulsar jKouveliotou et al.|1996^. 


One of the first properties of type II bursts discovered in the 
RB is the relaxation-oscillator behaviour: the fluence in a type II 


burst is proportional to the waiting time to the next burst (Lewin 
11976^. This is not observed in the BP ([Kouveliotou et aL 


et al. 


1996 1 . The relaxation-oscillator behaviour seems to point to the 


presence of a mass storage from which matter can be released and 
then replenished at a constant rate, so that the amount of depleted 
material determines the time to refill the buffer ( [Lewin et al.|1993| >. 
The accretion disc is the most obvious candidate for this buffer and 
all proposed models involve one (e.g., |Walker|1992l|Spruit & Taamj 
|1993| (. The various time scales at different positions in the disc can 
help to explain the range of observed type II burst time scales, while 
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a key role in acting as a gating mechanism might he played by a 


magnetic field (e.g., recent simulations by Kulkami & Romanova 


20(^ D^Angelo & Spruit|2010[ [Romanova et al.|2012[|D’ Angelo 


& Spruit|20T2). 

Notwithstanding that models involving some role for the NS 
magnetic field seem promising, the strength of the magnetic field is 
unfortunately not known in the RB due to the lack of pulsations. For 
the BP, estimates of the surface strength of the dipolar compo nent 
of the magnetic field B put it in the range (2 — 6) x 10^° G jDe-| 
Igenaar et al.|2014|l. While this is larger than in typical NS LMXBs 
(|Lewin et al.|19931l, it is not unique among bursters (estimates put 
the magnetic field of the NS in IGR J17480-2446 at 10® - 10“ G; 
e.g., ICavecchi et ^|2011^ . Therefore, at least another parameter 
besides the magnetic field strength (e.g., the stellar spin, or the in¬ 
clination of B with respect to the stellar spin axis) must determine 
the type II burst occurrence. 

[Bagnoli et al.| pOT^ hereafter B13) have shown that ther¬ 
monuclear (type I) X-ray bursts in the RB occur over an unusually 
large range of accretion luminosities, from about 45 per cent of the 
Eddington luminosity L^dd (assuming it is located in the globu¬ 
lar cluster Filler 1 at a distance of 7.9 ± 0.9 kpc, [Valenti, Ferraro| 
[& Origlia|2010[ l down to almost quiescence ( [Masetti et al.|2000[ ). 
However, while the type I bursts span the entire range of luminosi¬ 
ties, type II bursts only appear below a critical luminosity of about 
10 per cent Z/Edd ( [Guerriero et al.|1999[ B13). Hence, there are a 
type-I-only phase and a mixed phase. 

The type-I-only state is a soft state because the persistent 
emission spectrum is dominated by a disc black-body component, 
and it spans the range 10-45 per cent of Z/Edd- This is almost 
unique among type I burst sources: all other bursters show nearly 
no type I bursts above L ~ O.ll/Edd ( [Cornelisse et al.|2003^, ex- 
cept for the slow (11 Hz) pulsar IGR 117480-2446 l Linares et al.[ 
[2012[ >. The abundance of type I bursts in the RB above this thresh¬ 
old prompted B13 to propose that it is a slow rotator too, and to 
speculate that this might actually be one of the necessary ingredi¬ 
ents for producing type II bursts. 

As the RB exits the soft state, the first type II bursts appear. 
These first type II bursts are the most energetic, long and infre¬ 
quent. A brief (typically a few days) intermediate state exists in 
which these are accompanied by double-peaked type I bursts and a 
QPO at 0.25 Hz ( [Bagnoli et al. 12014] hereafter B14). The unusual 
shape of the type I bursts found in this short-lived intermediate out¬ 
burst phase seems to be due to a temporary obscuration of the NS 
surface during the burst decay, producing the double peaks. The 
ensuing type II bursts, the obscuration phenomenon and the occur¬ 
rence of a QPO are probably related to the change in the accretion 
geometry that is widely thought to take place between the soft and 
the hard state (e.g., review by [Done, Gierlinski & Kubota[2007[ and 
references therein). 

Having studied the soft and intermediate states in our afore¬ 
mentioned papers, we now complete our study of the entire Rossi X- 
ray Timing Explorer (RXTE, 1995-2012) RB database by focussing 
on the type-II-burst dominated hard state, so called because of the 
power-law component that dominates the persistent emission spec¬ 
trum. We aim at constructing a population study of type II bursts 
to compare it with the predictions of the few available theoreti¬ 
cal models, and as a future reference. In Sec. we introduce the 
dataset, and in Sec. we illustrate our methodology. In Sec. we 
present the type II burst properties and their relationship with the 
persistent emission and the type I bursts, taking advantage of the 
unprecedented detail provided by the RXTE data. We summarize 
our findings, some of which were previously unknown, in Sec. 


and discuss their physical implications in Sec.|^ Finally, in Sec.|^ 
we review proposed theoretical models, checking their predictions 
against the properties of our sample. 


2 OBSERVATIONS 

2.1 The general sample and burst identification 


This study is about measurements with the Proportional Counter 
Array (PCA) onboard RXTE. It consists of five co-aligned propor¬ 
tional counter units (PCUs) that combine to a total photon collect- 


ing area of 8000 cm over a 2 to 60 keV bandpass (Jahoda et al. 
|2006| (. 


Nearly all RXTE observations of the RB (exposure time 
2.4 Ms) were carried out in one of three configurations. For 553 ks, 
the PCA was directly pointed at the source. These observations are 
contaminated by the presence of the persistently bright 4U 1728-34 
at an angular distance of 0.56°. For 612 ks, the pointing was offset 
by just this amount, thus avoiding the contamination at the expense 
of roughly half the effective area. Finally, for 1.3 Ms the PCA was 
pointed at 4U 1728-34, serendipitously encompassing the RB. We 
have investigated all of the 2.4 Ms of data. 

We extracted a lightcurve from the ’STANDARD- 1’ data for 
all the observations containing the RB in the field of view (FOV). 
These data have 0.125 s time resolution and no photon energy res¬ 
olution. The telescope orientation was determined with the FTOOL 
PCACLRSP, and deemed stable if the jitter was below 20 arcsec. 

We performed a visual inspection of all the light curves. For 
the details of source (RB vs 4U 1728-34) and burst (type I vs type 
II) identification, we refer the reader to Appendix [A| In total, we 
identified 123 type I and 8458 type II bursts from the RB. 


2.2 Type II bursting modes 

As already mentioned, type II bursts are only emitted by the RB 
at persistent luminosities below 10 per cent of the Eddington limit 
(B13, [Guerriero et al.|1999^ , when the RB is in the hard state (B14). 
Actually the RB almost never reaches higher luminosities during 
outbursts after 1999 ( [Masetti|2002[ B14). 

Three modes of type II bursting behaviour follow one another 
in a smooth transition, differing in lightcurves, energetics and burst¬ 
ing patterns ( [Marshall et al.| 197^ [Guerriero et al.|1999^ . To illus¬ 
trate the evolution of the type II burst behaviour, we plot the light 
curve from RXTE data of the outburst in June-July 2006 (Fig.[^. 

The first type II bursts are hundreds of seconds long, often 
with flat peaks, and have the longest recurrence times (sometimes 
longer than a few ks observation). The persistent emission drops 
right before and after them. These are the so-called mode-0 type II 
bursts (Fig. ^ top-left inset). We identified a record long mode 0 
burst with a duration of at least 1100 s. It is visible in Fig.[T]as the 
long plateau on June 26. 

At lower persistent fluxes, the large bursts are increasingly 
preceded by sequences of 8 to 40 short bursts. After the final longer 
burst, a longer gap is measured, during which the persistent emis¬ 
sion rises to a level higher than that of the intra-burst emission, 
in a “hump” shape. This bursting pattern is called mode 1 (Fig.[^ 
top-right inset). 

As the outburst decays further, the sequences of mode 1 type 
II bursts become increasingly long, the final bursts less energetic, 
and the humps shorter and less pronounced. Eventually all bursts 
have similar energetics and durations: this is the so-called mode 2 
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outburst June-July 2006 



06/2CE2 23 24 25 26 27 28 29 30 07/01 02 03 04 05 06 07 08 09 10 11 12 13 14 15 


Figure 1. Sandwiched light curve of the observations of the outburst starting on June 20 2006, for the complete PCA bandpass at 1 s resolution resolution. 
Each consecutive data stretch is typically about 3000 s long. Data gaps mark jumps in time when there were no observations. Photon count rate per PCU is 
on the left y-axis, unabsorbed bolometric flux on the right y-axis. The horizontal dashed line indicates the Eddington luminosity, with the grey-shaded area 
indicating the uncertainties on the distance and the spectral fit (see Sec. ED- Count rates are not background subtracted in this plot, but are corrected for the 
number of active PCUs and the collimator response, (throughout the rest of the paper, the background counts have been subtracted). 4U 1728-34 was never in 
the field of view. The non-zero count rates before the onset of the outburst are therefore due to the instrumental background and the galactic ridge emission 
(see Fig. 2 in |Valinia & Marshall|1998} . The insets show the vaiying (following in time left to right, top to bottom) type II bursting modes during the outburst 
decay. Each inset shows 1 ks of data, plotted logarithmically between 50 and 5000 c s“^PCU“^. This outburst shows the prototypical sequence mode 0 
mode 1 mode 2 in bursting behaviour. The soft state, however, does not appear at the beginning of the outburst, hence preceding the type II bursting phase 
as in most cases, but it is briefly entered by the RB around June 27, giving rise to the long flat peak around that date. Note that the shorter plateau visible on 
June 26 is actually the longest type II burst in our sample, at about 1100 s (see Sec.|4.2). 


(Fig. □ bottom-middle inset). Finally, even these bursts become 
increasingly sparse, as the source moves to quiescence (Fig. 
bottom-right inset). 

We caution that the sequence “no type II bursts —>■ mode 0 
—>■ mode 1 —>■ mode 2” is only a general trend. For instance, the 
flat peak in Fig. [T] starting on June 27 is not a very long mode 0 
type II burst, but an actual switch of the source to the soft state, 
during which a type I burst can be seen, but no type II bursts. Two 
observations were taken the day before and the day after showing 
long mode 0 type II bursts, meaning that the soft state lasted briefly. 

Also, the RB can sometimes be seen switching back and forth 
between modes. The switch is sudden rather than gradual (as shown 
in Fig.|^ already noted by |Lewin et al.|1993) >, although the RB can 
be seen to transit back and forth between modes for a few ks. It is 
clear from the picture that the minima before and after a mode 0 
burst are not “dips”: rather, the humps in between are what is pe¬ 
culiar and are a feature common to all very energetic type II bursts, 
independent of bursting mode. 

The sampling along an outburst is often patchy, meaning that 


it is not easy to establish an average duration for each mode. The 
outburst in Fig.[^is probably the best sampled one. It features the 
longest mode 0 phase in the entire dataset, lasting ~ 8.7 days, albeit 
including a brief return into the soft state. The longest uninterrupted 
mode 0 bursting phase occurred in an earlier outburst, and lasted 
~ 6.5 days. Mode 1 bursts are visible in Fig.[T]for « 2.1 days. In a 
later outburst, they were observed to last for twice that time, « 5.2 
days. Finally, mode 2 burst were continuously observed for up to 
24 days, the longest stretch in the dataset. 

2.3 The type II burst sample 

Our search routine found 8458 type II bursts in the entire dataset 
(see Sec. ED- The results are summarized in Table [T] We deter¬ 
mined the instrumental and cosmic diffuse background with the 
FTOOL PCABACKEST. No bursts were analyzed that are outside so- 
called good time intervals (GTI), are affected by data gaps are over¬ 
lapping with type I bursts (from either the RB or 4U 1728-34), or 
are in ObsIDs for which the collimator response was not constant 
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4 Bagnoli et al. 



Figure 3. Histograms of all persistent flux parameters (see Sec. 4.l|: the persistent count rate immediately before a burst the average persistent count 

rate after a burst (Fpers) and the total average flux (including the burst) (Ftot)- The peak of the distribution, 5th and 95th percentiles are indicated for each 
parameter. Only bursts that were not contaminated by 4U 1728-34 have been included (see Sec. ED- 
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Table 1. The sample of analyzed^^^ type II bursts. 
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Figure 2. Examples of mode switching, 0^1 (top panel) and 1^0 (bot¬ 
tom panel). Time axis with respect to the start of the ObsID. The y-axes 
are as in Fig. The horizontal dashed line indicates the Eddington lumi¬ 
nosity, with the grey-shaded area i ndic ating the uncertainties on the dis¬ 
tance and the spectral fit (see Sec. 3.l|. The dotted line shows = 

60 c s“^PCU“^. The quiescent level between quick bursts is the same 
as the dip level before and after mode 0 bursts. The RB therefore does not 
“dip” but rather shows humps of enhanced emission after long type II bursts. 


because of a slew. This left us with 7601 burstsj^Of these, 3662 
type II bursts are present in 389 ks of direct pointings, while 3183 
were observed during 341 ks of offset pointings, and 756 during 
92 ks of 4U 1728-34 observations. 

In total, the RB was stably in mode 0 during the entirety of 
an ObsID for 114 ks, in mode 1 for 170 ks, and in mode 2 for 
492 ks. In 51 ks of data, the behaviour was intermediate between 
two consecutive modes, switching often back and forth during the 
transition. We found 104 mode-0 bursts, 2522 mode-1 bursts and 
4971 mode-2 bursts. For 4 bursts taking place during very short 
observations (a few hundred seconds) we could not clearly assign 
a mode. 


^ A table including all the burst properties for the entire sample is available 
online. 


Bursting mode^^^ 

Number of bursts 

RB offset 

per aimpoint^^^ 
1728-34 

Total 

0 

78 

24 

2 

104 

1 

964 

895 

663 

2522 

2 

2620 

2260 

91 

4971 

unclear 

0 

4 

0 

4 

alt 

3662 

3183 

756 

7601 


About 800 bursts identified by our routine were not analyzed, see 
Sec.O 

(b)For a description of the bursting patterns, see Sec |2.2| 

^^^The different orientations of the telescope in the data and their effect on 
FOV contamination and effective area ai‘e discussed in Sec. o 

Furthermore, during 210 ks of data, the RB was at such a high 
luminosity that only type I bursts were observed (B13). Finally, no 
RB bursts of either type were found in ObsIDs totalling 1.4 Ms 
(mostly in observations with the instrument pointed at 4U 1728- 
34), when the RB was in quiescence. 


3 METHODOLOGY 

3.1 Counts as proxy for energy 

We adopted the photon counts as a proxy for the energy flux, rather 
than fitting time-resolved spectra. This approach offers an obvious 
advantage in speed of the data reduction. Furthermore, the intra¬ 
burst persistent emission would often be too short and weak to 
constrain the spectrum well. Finally, the flux variations, which can 
be very large during type II bursts, are largely achromatic; it is 
mostly the normalization factor and not the intrinsic spectral shape 
that is changing. This is shown in Appendix]^ where we perform 
time-resolved spectroscopy on a representative selection of the type 
II bursts. We obtained the best fits with a Comptonization model 
(COMPTT in XSpec |Titarchuq|1994| >. The conversion factor be¬ 
tween PCA counts and 3-25 keV fluence is 1.02x 10“^^ erg cm“^ 
per c/PCU. The uncertainty induced by the spectral shape is 10 
per cent (full range, not la). The bolometric correction would be 
1.24 ± 0.07 for the COMPTT model. 

For a “canonical” NS (mass M = IAMq) the Eddington lu¬ 
minosity Lfidd is 2.lx 10®® erg s“®. This assumes isotropic emis¬ 
sion, no gravitational redshift and an accreted H fraction X = 0.7 
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because the RB type I bursts show long (100 s) tails that are typical 
of burning of H-rich material. For a distance of (7.9 ± 0.9) kpc, 
this yields a flux FEdd = (2.8 ± 0.3)xl0“® erg cm“^ s“^. 
Given the above conversion to the bolometric flux, and propagat¬ 
ing the relative errors on the spectral shape, bolometric correction 
and distance, the Eddington luminosity in the RB corresponds to 
Fidd = (2200 ± 500) c s-^PCU"^ 

3.2 Extracting the sample properties 

For each type II burst, we calculated from the STANDARD-1 light 
curve the start time and duration fdur at 0.125 s accuracy by seek¬ 
ing backward and forward in time for the first 2 consecutive time 
bins where the count rate was below the persistent emission level. 
This way, we also determined and At"®’'*, the recurrence 

times with respect to the previous and the next burst respectively, 
as the difference between the start times. Because the burst dura¬ 
tion can be a significant fraction of the duty cycle for type II bursts 
(unlike for most type I bursts), the actual non-bursting time can be 
significantly shorter than the recurrence times so defined. 

We defined two different measurements of the intra-burst per¬ 
sistent flux. First, we calculated the actual intra-burst fluence Fpers 
and divided that by the intra-burst time, providing the average per¬ 
sistent flux (Fpera). lu a vast majority of cases the instantaneous 
flux was stable during the interval, and deviated little from the av¬ 
erage. 

As mentioned in Sec. |2.2| the persistent emission however 
tends to rise after very energetic bursts (the mode 0 and end-of- 
sequence mode 1 bursts), reaching a minimum immediately be¬ 
fore and after those bursts. Therefore, for the purpose of measuring 
the net burst emission properties, we also calculated the average 
count rate in the 10 s before a burst (or less, when AF'^®'' is 
shorter). This assumes that, during the burst, the persistent emission 
will be stable at the minimum level reached right before it. 

We made this choice for two reasons. Firstly, this assumption 
has long been part of the standard approach to burst analysis (e.g., 
IGalloway et al.|2008^ . Secondly, the reappearance of the persistent 
emission at roughly the same flux after the energetic bursts (see 
e.g.. Fig. supports our choice. 

We subtracted from the observed count rates to calculate 
the net burst emission properties. We defined the peak flux Fpeak 
and the total burst fluence Fb as the maximum net count rate and 
the integrated burst net counts during this interval. 

The ratio of Fpers to Fb is the parameter a, which is key to 
distinguishing between type I and II bursts, because the ratio of the 
fluence liberated in the persistent emission to that the burst diag¬ 
noses whether the latter can be explained by thermonuclear burning 
of the accreted material or not. For reasons that will become clearer 
in Sec. |4.2[ we defined fwo values of a for each bursf, employ¬ 
ing Fpers before and after fhe burst: thus we retrieved a backward- 
defined qP’^®'' and a forward-defined a"®**. The previous is the pa¬ 
rameter that is normally employed in studies of type I bursts, but 
as we will see, the latter shows a more significant relationship for 
type II bursts. 

We subtracted the instrumental background counts (calculated 
and interpolated from Standard-2 data in 16 s resolution) from 
all quantities, and then normalized them by the number of active 
PCUs and the correction given by the telescope orientation (see 
Sec. |2. l| l. Clearly, all quantities related to the persistent emission 
have to be considered upper limits (lower for a) when the FOV 
is contaminated by 4U 1728-34. We did not include these in our 
analysis. 



Figure 4. The average values of (Fpers) and (Ftot) during mode 

0, 1 and 2 type bursts (see Sec. |4. 1) . The error bar indicates the standard 
deviation in each sample. 


4 RESULTS 

4.1 The persistent emission 

First we look at the properties of the persistent emission when the 
RB is emitting type II bursts, whenever not contaminated by 4U 
1728-34 (see Sec. \2.1\ . We plot the various properties in Fig. 
with the peak of the distribution (the highest bin) and the 5 and 95 
percentiles indicated. 

The average persistent emission count rate before a burst 
is limited to a narrow range of values, with 90 per cent of the distri¬ 
bution between 25 and 60 c s“^PCU~^, or ~ (0.01 — 0.025)FBdd. 
Therefore, no bursts are observed on top of a null net persistent 
emission. The largest observed value is 110 c s“'^PCU“'^, corre¬ 
sponding to ~ 0.05FBdd (note that at times when there are no type 
II bursts, the persistent flux can be as high as ~ 0.45FEdd). 

As explained in Sec. |3.2[ this is only a lower limit on the 
actual intra-burst emission. Flowever, short mode-1 and mode-2 
bursts make up the majority of our sample, so that the follow¬ 
ing histogram, of the average persistent count rates after a burst 
(Fpers), shows a very similar distribution, with a slightly longer 
tail. The minimum values are similar to those of F^J." because 
the persistent emission never dips between bursts, but only right 
before and after them. The maximum of the distribution is larger 
with respect to F^)", with all but ten bursts showing (Fpers) < 
200 c s-^PCU“^ ~ 0.09FEdd. 

Since both the persistent emission and the type II bursts prob¬ 
ably originate from the release of gravitational energy (although 
possibly through different channels), the total average flux in the 
interval between two bursts At, including both the burst and per¬ 
sistent emission, is also a potential proxy for the actual mass accre¬ 
tion rate of the RB during the hard state (under the assumption that 
the radiative efficiency is similar for the persistent emission and the 
type II bursts). We call this flux (Ftot), and plot it in the following 
histogram in Fig.[^ The bursts seem to increase the accretion flux 
by a factor of ~4, up to an upper limit of 400 c s“^PCU“^, or 
~ 0.17FEdd. Therefore, an average 75 per cent of the accretion 
goes through bursts. 

All three measurements show a monotonically decreasing 
trend from mode 0 to 1 to 2, as plotted in Fig. Flowever, the 
dispersion within each sample is similar or larger than the differ¬ 
ence in values between modes. We obtained the same result try- 
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Figure 5. Histograms of all burst parameters. The peak of the distribution, 5th and 95th percentiles are indicated for each parameter. The dashed line in the 
left-top plot is the flux corresponding to the Eddington luminosity for M = 1.4Mq, X = 0.73. The grey area reflects the uncertainty in the spectral fit (see 
Appendix [B). The shortest observed Af = 2 s, = 0.125 s. The dashed line in the a plots divides the regions where the burst fluence exceeds that of the 
persistent emission and viceversa. ol values have only been computed for observations that were not contaminated by 4U 1728-34 (see Sec. lin¬ 


ing to reproduce the same figure for the outburst plotted in Fig. m 
alone. It seems therefore that the persistent emission variations only 
marginally trace the changing type II burst behavior. 


4.2 The type II bursts 

For all 7601 analyzed type II bursts we study their energetics (peak 
flux Fpeak and fluence Eh) and timescales (duration tdui and re¬ 
currence time At). Furthermore, we investigate their relationships 
to the persistent emission properties for the 3183 bursts observed 
when 4U 1728-34 is outside the FOV. In Fig. we have plotted 
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histograms of the burst properties, again annotated with the peak of 
the distribution (the highest bin) and the 5 and 95 percentiles. 


4.2.1 Burst properties 

Starting with the distribution of peak fluxes of type II 
bursts, the peak is at about half the Eddington luminosity, or 
1100 c s“^PCU~^. Overall, the distribution is relatively nar¬ 
row, with 90 per cent of the bursts peaking in the range (1.1 — 
2.7)X10® erg cm“^ s”'^. The distribution seems compatible with 
being Eddington limited. Only 1 per cent of the bursts are above 
the upper limit on Z/Bdd given by the uncertainty over the distance 
and the spectral fit (see Sec. [m this limit is actually more clear in 

Fig-IZig. 

A larger dynamic range of values is observed for the burst 
fluences, which span about three orders of magnitude, making for 
a fairly large distribution. The majority of bursts are in mode 1 
or 2, so it is not surprising that many show a relatively small 
fluence. These bursts radiate typically between 10“® and several 
10“^ erg cm“^. The more energetic mode 0 bursts show mini¬ 
mum fluences of several 10“^ erg cm“^. 42 mode 0 bursts show 
i?b >10“® erg cm“^, peaking at = 5.8x10“® ergcm“^. 

The durations of type II bursts can vary by 4 orders of 
magnitude, although 90 per cent of them last between 2.4 and 
23.75 s. The shortest bursts actually challenge the time resolution 
of Standard- 1 data of 0.125 s, and need to be resolved in event¬ 
mode data (Eig.|^ top). The shortest burst in our sample only lasted 
130 ms. We found 61 bursts with sub-second durations, all during 
mode 2 bursting phases towards the end of an outburst, 46 of them 
in ObsID 92026-01-08-05 alone. 

These very short bursts do not show peak fluxes smaller than 
the rest of the sample. Their fluences are however some of the 
smallest observed, liberating on average ~ 5x10“® erg cm“®. 
Their light curves tend to show very fast rises and equally fast de¬ 
cays, although they can have shapes that are unseen in the rest of 
the sample (such as the burst in the bottom-left panel of Fig. 
which has a gradual rise and a sudden decay). 

The longest burst observed in its entirety lasts 533 s (Fig.|^ 
bottom). The record holder is however a burst lasting at least 
1100 s, the rise of which was not caught in the observation. Only 
26 bursts recorded durations longer than 100 s, all but two being 
mode-0 bursts. 

The recurrence times are mostly below 200 s (again, with 
the exception of the mode-0 bursts), the shortest and the longest 
(without data gaps) At being 2 s and 1656 s respectively. Over¬ 
all, the distribution is extremely asymmetric. This is logical, since 
the shorter the recurrence time, the more bursts will be observed. 
The peak of the distribution is however well above zero, at At — 
15 — 18 s. Below this bin, the distribution drops abruptly. Concern¬ 
ing mode 0 bursts, their recurrence times (excluding lower limits) 
range between 105 and 1562 s. 

For the 3183 bursts that were in observations uncontaminated 
by 4U 1728-34 we compared the fluence released in the burst to that 
in the persistent emission. We calculated the and val¬ 

ues, respectively, by dividing the integrated fluence in the interval 
prior and after the burst by the total burst fluence Eh. The former 
is the parameter that is employed when dealing with type I bursts, 
since thermonuclear bursts are due to the burning of material after 
this has been accreted. Fig. shows that the values of for 

type II bursts peak at 0.20, meaning more energy is released in the 
bursts than in the persistent emission. They are incompatible with 
the values expected from thermonuclear burning (~ 40 at least, for 
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Figure 6. (Top.) Examples of some of the shortest bursts in the sample 
from the 10 ms resolution light curve of ObsID 92026-01-08-05 (July 20 
2006). Left y-axis is as in Fig.^ The bursts take place during an otherwise 
ordinary sequence of mode 2 type II bursts. (Bottom.) For comparison, the 
longest entirely observed burst in the sample, at 533.125 s, from ObsID 
40058-02-03-0IS (April 20 1999). The y-axes are as in Fig. The flat 
top is substantially below the Eddington limit, indicated by the horizontal 
dashed line, with the grey-shaded area indicating the uncertainties on the 
distance and the spectral fit (see Sec. HD- A longer burst (> 1100 s, see 
Fig.[T]and Sec.|4.2j was only partially observed. 


H rich fuel, |Lewin et al.|1993| and the minimum measured a ~ 10, 
IKeek et al.|20T0) . 

Even at the lowest persistent fluxes, when the RB approaches 
quiescence and the type II burst frequency drops, (see bottom- 
right inset of Fig. 10. The only type II bursts with 

larger qP®®'' values are those with sub-s durations, which can have 
aP®®'' ~ 20 — 200, but are much shorter than type I bursts. 

The type II bursts are thought to draw from the same en¬ 
ergy reservoir as the persistent emission. Also, their occurrence 
and properties seems to affect the persistent emission that follows, 
rather than to depend on the persistent emission that precedes them 
(see below). Because of this, we also plot the integrated ratio of the 
fluence/o/Zowing a burst to the integrated burst fluence, a®®’'*. The 
two distributions are similar, although that of a“®’‘* is narrower and 
does not approach values as small as the other one. The final large 
bursts in mode 1 sequences tend to have short Afprev (the time in¬ 
terval to the previous burst), and long Afnext. Hence their qP'^®'' 
are smaller than their o?"®’'*. 
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Figure 7. Peak flux .Fpeak burst duration a function of the burst 

fluence E^, (see Sec. |4.2^ . Count rates and energies are as in Fig. ^ The 
horizontal dashed line indicates the Eddington luminosity, with the grey- 
shaded area indicating the uncertainties on the distance and the spectral fit 
(see Sec. in}- The correlation (given by the Spearman correlation factor p) 
is strong in both cases, although it saturates in the upper plot at the flux 
con'esponding to the Eddington luminosity. The dashed-dotted line in the 
bottom plot indicates the duration of an Eddington-limited burst of a given 
fluence. The red diamonds are the values predicted by three the models in 
Fig.[T5] with two of them nearly overlapping (see Sec. ED- 


4.2.2 Correlations 

Having described the burst properties individually, we turn our at¬ 
tention to the relationships between them. As we show in Fig. |7] 
some burst properties show trivial correlations. The burst fluence 
Eh correlates with both the peak flux Fpeak (Spearman rank-order 
correlation coefficient p = 0.57) and, more strongly, with the burst 
duration fdur (p = 0.90), with the mode 0 bursts having the largest 
Eh and longest tdui. In the lower left comer of the upper plot 
are the bursts for which almost the entire fluence is emitted in a 
short interval around the burst peak (in the lower plot these bursts 
are mostly shorter than 10 s). At the opposite end of the top plot, 
the correlation between Fpeak and Eh saturates once the peak flux 
reaches the Eddington limited. There seems to be no correlation 
between peak flux and duration (not plotted). 

It is also interesting to see how these quantities, which deter¬ 
mine the burst energetics, depend on the bursting rate at a given mo¬ 
ment. Fig.j^shows the type II burst peak flux, fluence and duration 
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Figure 8. Burst peak flux Fpeak^ fluence E'b’ duration f^ur us a function 
of the waiting time At to the next burst (see Sec. |4.2) . Count rates and 
energies are as in Fig.^ The horizontal dashed line indicates the Eddington 
luminosity, with the grey-shaded area indicating the uncertainties on the 
distance and the spectral fit (see Sec. ED- The middle and bottom plot show 
the relaxation oscillator behaviour of the type II bursts; the burst recurrence 
time seems to more accurately predicted by the duration of the previous 
burst than its fluence. The least-squares fit to the data in the bottom plot is 
also shown with a dotted line (see Sec.|4.2^. 
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Figure 9. Same burst properties as in Fig.j^ here plotted against the average 
persistent flux (Fpers) (see Sec. |4.2| for definition). Only bursts observa- 
tions that were not contaminated by 4U 1728-34 are plotted (see Sec. HD- 
Count rates and energies are as in Fig. The horizontal dashed line in¬ 
dicates the Eddington luminosity, with the grey-shaded area indicating the 
uncertainties on the distance and the spectral fit (see Sec. ED No overall 
con'elation is visible in the data, although it might hold for the longest, most 
energetic bursts (see Sec.|4.2). 


[10’’^ erg/cm^] 



E[;ump = Epers-(AtnextXFp™?3) [c/PCU] 


Figure 10. The burst fluence E^, against the “net” (i.e. above hump 

fluence Bhump (see Sec. |4.2[ . Count rates and energies are as in Fig. 
Above a critical ^ 2 X 10"^ c/PCU, a correlation appears. Both the 
brightest mode 1 bursts and all mode 0 bursts can be above this threshold. 


as a function of the recurrence time to the next burst On 

top, one can see that Epeak is roughly independent of the waiting 
time (p = 0.04). A correlation instead appears in the middle panel, 
which plots Eh against At"®’'* (p = 0.71): this is the well-known 
relaxation-oscillator behaviour of the RB, according to which the 
more energetic a burst, the longer it takes for the next one to occur. 
This is the reason why we chose to introduce the parameter q"®’^* . 
However, a look at the bottom plot in Fig. [^reveals that Af“®’‘* 
correlates more tightly with fdur (p = 0.87) than with Eh- This 
implies that the duration of a type II burst, rather then the energy 
it releases, determine the time it takes for the next type II burst to 
occur. A least-squares fit yielded Fee oc There 

is, however, a large spread in the data, especially for the mode 1 
and mode 2 bursts. 

As we explained in Sec. |2.2| different patterns of type II bursts 
follow each other during the outburst decay, while on average the 
persistent flux decreases. We therefore wish to assess whether the 
burst properties show any dependence on the strength of the persis¬ 
tent emission. Fig.|^shows the burst peak flux, fluence and duration 
plotted against the average persistent flux (Epers)- It is immediately 
evident that all correlations are weaker than with respect to Af®®’‘‘. 
As in the previous case, Epeak shows no evident relation to the per¬ 
sistent emission properties (p = —0.03). In the middle and bottom 
panels, where we plot Eh and fdur against (Epeia), no monotonic 
trends can be seen. Instead, two distinct behaviours appear. While 
for bursts that liberate little fluence (Eh < 2 x 10"^ c/PCU) or that 
are short (fdur 20 s) Eh and fdur do not correlate with (Epers) 
(p = —0.04 in both cases), for energetic and long bursts a positive 
correlation is visible, that at least for Eh appears to be significant 
(p = 0.46 and 0.23 respectively). 

In other words, when looking at how the burst properties are 
influenced by the persistent emission, there is a dichotomy in be¬ 
haviour between the mode-0 type II bursts, which have large flu- 
ences and long durations, and appear in the upper half of the plots, 
and mode 2 bursts, characterized by small fluences and short dura¬ 
tions, appearing in the bottom half. Mode 1 bursts can clearly be 
found both in the lower half of the plots of Fig. (low fluence and 
short duration) and in the upper half (large fluence and long dura¬ 
tion). These are, respectively, the short bursts initiating the mode 
1 sequence, and the long one that concludes it. These two groups 
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Figure 11. Light curve and hardness ratio of ObsID 50420-01-09-00R, fea¬ 
turing a type I burst taking place during a mode 0 type II burst. The vertical 
lines indicate the approximate start and end time of the burst. The net peak 
flux to persistent flux ratio ^ < 0.2, hiding the hardening of the flux during 
the type I burst below the detection level jLinares et al.|20lT] B13). 


I 

i 

I 



Figure 12. Light curves of 5 type I bursts in the hard state plotted against the 
average profile of the 22 single-peaked type I bursts featuring the longest 
and brightest decays (see B14), all of which are instead observed in the soft 
state. The net bursts counts have been normalized so that the average flux 
during the burst peak (which is defined as the interval during which the flux 
is above 90 per cent of the burst peak flux) is 1. The vertical extent of the 
shaded region indicates the standard deviation of the soft-state sample. The 
hard-state type I bursts were randomly chosen among those taking place 
on top of a stable persistent emission, i.e. outside humps, and showing no 
overlapping type II bursts during the first 40 s of the type I burst. 


inhabit the same regions of the plots as the mode-2 and 0 bursts 
respectively. 

The same threshold in fluence seems to determine the presence 
of a hump after a burst. Fig. [T^ shows that, independent of burst 
mode, above a critical E\, ~ 2x 10“^ erg cm“^, the burst fluence 
correlates with the “net” (i.e. above hump fluence JShump- 

Again, bursts above this threshold are not only the mode 0, but also 
the resembling mode 1 that end a sequence of bursts and precede a 
hump. 


4.3 Type I bursts in the mixed (hard) state 

Various methods help identify the type I bursts in the sample. Their 
profiles in the light curves differ for a number of features from those 
of type II bursts. The type-I burst peak fluxes are smaller than those 
of most type II bursts (type I bursts in the RB reach at most one third 
of the Eddington luminosity, B13). Also, all type II bursts show 
multiple peaks in their decays, while type I bursts generally have 
smooth decays (but see B14). Finally, type I bursts have generally 
longer durations, causing the fast-recurring type II bursts to overlap 
with their exponential decays. Finding candidate type I bursts is 
therefore relatively straightforward. 

Nonetheless, a secure identification of type I bursts should 
come from measurements of a or from the detection of cooling 
in their time-resolved spectra. However, the first method is only 
applicable in a state where either type I or II bursts are emitted. 
In a mixed state, where type II bursts recur faster than the type I 
bursts, one cannot rely on their Oprev values to discriminate be¬ 
tween them. The second method, the detection of cooling, is only 
a sufficient, and not necessary condition for the type I burst iden¬ 
tification. Observing spectral cooling is generally only possible for 
bright type I bursts (Fpeak ^ 0.Topers; [Linares et al. p0TT] >. B13 
also showed that some RB bursts which lack evidence for cooling 
are to be identified as type I, due to their large a values and their 
occurrence in the brightest phases of the soft state. 

To determine the presence of cooling, we measured hardness 
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Figure 13. Two examples of type I bursts during the hard state, during mode 
1 (top) and 2 {bottom) type II bursts respectively. The x axis shows time 
with respect to the type-I-burst start time tQ. The y-axes are as in Fig.^ The 
horizontal dashed line indicates the Eddington luminosity, with the grey- 
shaded area indicating the uncertainties on the distance and the spectral fit 
(see Sec. ED Notice that the type I burst during the mode 1 phase takes 
place during the hiatus in type II bursts that follows a long one, while the 
persistent emission level rises in a hump shape. This is the case for 10 out 
of 12 type I burst during the mode 1 phase Neither effect is observable as a 
type I burst takes place in the mode 2 phase (see Sec. |4.3) . The type I burst 
during the mode 2 phase briefly affects the type II burst emission. This is 
only obseiwed in 1 case out of 11. 
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ratios as the flux ratios in the 5.1-7.6 keV and 2.3-5.1 keV bands. 
26 candidate bursts show secure evidence of cooling, and we can 
therefore confirm they are type 1 even without being able to mea¬ 
sure their a values. 

Additionally, we found three bursts taking place “on top” of 
a mode 0 type 11 bursts (see Fig. [TTJ. For these bursts Fpeak < 
0.2F'pers, meaning the lack of cooling in their spectra cannot rule 
out their type 1 nature. Given their distinctive shape (fast rise and 
exponential decay), these three bursts are unlikely to be a sudden 
spike of the type 11 bursts, which tend to show nearly flat plateaus 
(see Sec. |2.2^ . Therefore we believe they are type 1 bursts as well. 

Summarizing, we found a total of 29 type 1 bursts occurring in 
an observation where also type 11 bursts appeared, out of a total of 
123. Note that 17 candidate type 1 bursts could not be confirmed. 

As shown in Fig. [T^ the type 1 bursts in the hard state are 
significantly different from those in the soft state. Most remarkably, 
the durations are much shorter, with fdur ^40 s, as are the peaks, 
that lack the 6-8 s plateau displayed by type 1 bursts in the soft state 
(B13). 

Out of 29 type 1 bursts identified in the hard state, 6 took place 
around mode 0 type 11 bursts, 12 during mode 1, and 11 during 
mode 2. Ten out of 12 type 1 bursts found during mode 1 activity oc¬ 
cur during the intervals of enhanced persistent emission (“humps”) 
that are observed between two more tightly packed sequences of 
type 11 bursts (see Fig[^ top). Only two type I bursts are instead 
observed “inside” a sequence of mode 1 type 11 bursts. No such 
preferential location is observed for the 11 type 1 bursts occurring 
during the mode 2 phase, when the persistent emission count rate 
and the recurrence time of the type 1 bursts are roughly constant. 
However, it must be noted that one such type 1 burst seems to briefly 
(< 50 s) affect the type 11 emission (see Fig[T^ bottom), which is 
otherwise stable in this observation, as is common in mode 2. 

For mode-1 type-11 bursts, one can calculate the chance of ran¬ 
dom coincidence of a hump and a type 1 burst. We examined all the 
mode-1-burst sequences during the outburst plotted in Fig.^ and 
summed the integrated fluence from the beginning of each burst in 
a sequence until the last, longer one. We divided this integrated en¬ 
ergy by the fluence emitted during the following hump. Over the 
course of the outburst, this ratio decreases from 13.5 to 2.7. In or¬ 
der to reach a conservative estimate, we take the latter as the ratio 
of the amount of mass accreted during a mode-l-burst sequence to 
that accreted during a hump. This means that, for each individual 
burst, the probability of random occurrence of a burst on a hump 
is Psingie = 0.27. Au estimate of the probability that 10 out of 12 
type 1 bursts randomly take place on a hump rather than overlap 
with mode 1 bursts is thus given by the binomial distribution as 
Pio/12 = 7.24 X 10“®. Clearly type 1 bursts have a statistically 
significant preference for occurring during breaks in the type II ac¬ 
tivity, rather than overlapping them. 

B13 reported type I burst recurrence times in the soft state in 
the range 0.4-4 ks. As the RB enters the hard state, free becomes 
difficult to measure, as it typically surpasses the duration of an ob¬ 
servation. The numbers should therefore only be treated as upper 
limits. 

Nonetheless, we can establish average recurrence times of 
type I bursts across the different modes of type II bursting. There 
are 6 confidently identified type I bursts in 114 ks of observations of 
mode 0, 12 in 170 ks of mode 1, and 11 in 492 ks of mode 2. If one 
includes the dubious cases, the type I bursts are 13 in mode 0, 21 in 
mode 1 and 12 in mode 2. This gives upper and lower limits, respec¬ 
tively, on the average recurrence times: ~ (8.3 —18.1) ks, 

~ (8.1-14.1) ks and ~ (41.0-44.7) ks. Thus, 
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the type I burst recurrence time is roughly the same during mode 0 
and mode 1. 

We plot these recurrence times in Fig. [T^ together with the 
type I burst recurrence times measured in the soft state, against the 
persistent flux Fpers (B13). The latter was taken to be a proxy of the 
mass accretion rate on the source, and the data show a correlation, 


only measures the emission taking place between type II bursts, but 
these should also accrete nuclear fuel onto the surface. Therefore, 
we plot the type I burst recurrence times across the three modes 
versus both (Fpera) and ftot, which includes the burst fluence in 
the average. 

As Fig.[T^shows, regardless of the proxy that is chosen for the 
mass accretion rate, the largest recurrence time for type I bursts is 
observed during mode 2 type II bursts. Interestingly, observations 
line up with the relation extrapolated from the soft state when em¬ 
ploying (Fpers), while for Ftot there appears to be a scarcity of 
type I bursts in the hard state, across all modes. 


with free OC 


T The (Fpers) parameter we define in Sec. 


3.2 


5 OBSERVATIONAL SUMMARY 

The large sample of data we have studied allows us to address fun¬ 
damental questions over the type II bursts in the RB. In this section 
we first summarize the overall properties of the sample, starting 
with those that were already previously known: 

1) There are four patterns of type II bursts (e.g., Fig.[TJ. 

i) a (type II) burst-less soft state that is visible when the per¬ 

sistent emission flux is above 10 per cent of the Eddington 
limit. This occurs mostly at the beginning of outbursts. 

ii) During the outburst decay, the first type II bursts to appear 
are of so-called mode 0 ( |Guerriero et al.|1999^ . They have 
the longest durations and recurrence times, and the largest 
fluences. The flat tops often observed in long mode 0 bursts 
are often sub-Eddington. Indeed, the flatness of the bursts 
has before been noticed to be a property of long bursts jTan| 
|et al.|199T) , regardless of their luminosity. 

iii) The mode 1 burst pattern follows, with sequences of fre¬ 
quent short bursts ending with a longer and brighter one that 
is similar in appearance to the shortest mode 0 bursts. 

iv) Finally, the mode 2 pattern appears, consisting of short 
bursts like the short ones in mode 1 (but not interrupted by 
larger bursts followed by longer recurrence times). 

This is only a general rule and exceptions exist. For instance, 
the RB sometimes briefly switches back to the previous mode, 
and even back into the soft state. Mode 0 and 1 bursts can last 
for up to ~ 5 — 6 days, while mode 2 bursts persist for up to 24 
days. 

2) The persistent emission decreases monotonically as the RB pro¬ 
gresses from mode 0 to mode 1 to mode 2 bursts (Fig.[^. No 
type II bursts are observed on top of undetectable persistent 
emission. 

3) All bursts roughly follow a relaxation oscillator relation (Fig.[^ 
middle), meaning that the fluence in a burst determines the 
waiting time to the next burst. 

4) 95 per cent of the bursts show a values below 1.5, one order 
of magnitude below the minimum value that is possible for full 
thermonuclear burning. The distribution peaks at 0.2, meaning 
that 5 times more fluence is liberated in the bursts than in the 
persistent emission. This underlines the nature of type II bursts 
as accretion processes. 
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5) We calculated the probability of random occurrence of the ob¬ 
served fraction of type I bursts on a hump to be Pio /12 = 
7.24 X 10“®. Clearly, type I bursts during mode 1 have a sta¬ 
tistically significant preference for occurring during breaks in 
the type II activity, rather than overlapping them (see Fig|13|l. 
This preference was first noticed in SAS-3 data by| Ulmer et al.| 
\\911) , although it was [Hoffman et al.| jl978^ who first pro¬ 
posed that the bursts on the humps and those surrounding them 
were type I and II respectively. They also noticed that these 
bursts were “anomalous” with respect to the type I bursts taking 
place during what we call mode 2, in that they had smaller peak 
fluxes and sometimes lacked spectral softening during burst de¬ 
cay. Likewise, we could confirm the presence of spectral soft¬ 
ening in 11 out of 12 candidate type I bursts during mode 2, but 
only in 12 out of 21 during mode 1. 

New findings in RB type II bursts are: 

6) The range of minima in the persistent emission between bursts 
F^eSi is between 1 and 2.5 per cent of the Eddington luminos¬ 
ity. The distribution of the actual average persistent fluxes (i.e. 
taking the humps into account) (Fpers) extends to only slightly 
larger values, because the general sample is overwhelmingly 
dominated by the short-mode-1 and mode 2 bursts, without 
humps in between. Taking into account the burst fluence as 
well, the RB is on average accreting at up to ~ 0.17LEdd 
(Fig.[3|. 

7) The persistent flux does not generally seem to correlate with 
the burst properties (Fig.|^. Each burst mode presents a large 
dispersion in the range of observed persistent fluxes (Eig. |^, 
making the latter a poor predictor of the burst mode. 

8) The most energetic type II bursts reach the Eddington limit 
(^Edd = 2.8x 10“® erg cm“^ s~^; Eig.[^ top-left). 1 per cent 
of the bursts are super Eddington. However, since the maximum 
F^JSl = 3.0x10“® ergcm“^ s“^, even a slight decrease in the 
H abundance or moderate beaming are enough to bring every 
single burst within the uncertainty on I/Edd (Sec. [nj- Bursts 
reaching the Eddington limit are not accompanied by evidence 
for photospheric radius expansion in the time-resolved spectra. 
Generally speaking, the peak flux shows no correlation with 
the waiting time or the persistent flux properties, and the dis¬ 
tribution of peak fluxes is relatively narrow. Even the only visi¬ 
ble correlation, the one with the burst fluence (Eig.|^, shows a 
large spread. This suggests that the speed at which energy is re¬ 
leased does not depend on the total amount of energy available 
for a burst, nor on the global mass accretion rate. 

9) There is a clear dichotomy between mode 0 and mode 2 bursts. 
The former appear to feed back into the persistent emission, 
which rises in between them, while the latter do not seem to 
affect it. This is also reflected in the fact that only for mode 0 
bursts there appear to be correlations between some properties 
of the bursts and of the persistent emission (Figs. EE- Mode 1 
bursts clearly are an intermediate state of the instability, repro¬ 
ducing both behaviours. 

10) Long bursts - meaning both mode 0 bursts and the mode 
1 bursts that “close” a sequence of bursts - are followed by 
humps, i.e. burst-free enhanced persistent emission. This is par¬ 
tially contrary to what was previously believed, that the persis¬ 
tent emission “dips” between mode 0 bursts (Marshall et al.j 


fT^lKumeda et al.|1984||Stella et al.|7^|Lubin et al.|19921 

[Le 


etal.|1993[[Gu 


et al.|199^, despite the fact that 


there is no evidence for enhanced absorption in the spectra 
( [Stella et al.||1988[ [Lubin et al.[[1993^ . However, looking at 
periods when the RB switches between these modes of emis¬ 


sion, we clearly show that such “dips” are actually at the same 
minimum level observed during a sequence of mode 1 bursts 
(Fig.§. Since the persistent emission level is stable inside a 
mode 1 sequence, and no increase in the column density Nu 
is observed in the spectra of the “dips” (see Appendix]^, we 
argue that in both cases the emission following a long burst is 
enhanced, to decay again in the time running to the next burst. 

11) The recurrence time distribution of type II bursts is strongly 
peaked at At = 15—18 s, and drops abruptly below this 
value, while showing a smoother decrease to larger recurrence 
times (Fig. |^. Bursts with shorter recurrence times are found 
(as short as 2 s) but they are a much smaller fraction of the sam¬ 
ple, despite the fact that for a shorter At more bursts should be 
produced. This implies that a minimum time exists for the in¬ 
stability to develop, and that it is difficult for the RB to produce 
type II bursts more quickly than that. 

12) We found 61 bursts with durations shorter than 1 s, and as 
short as 130 ms (Fig|^ top), one order of magnitude smaller 
than the shortest bursts so far reported for type II bursts (2 s 
in EXOSAT data, [Lubin et al.||199iy This puts a conshaint 
on theoretical models for this instability, that appears to cover 
four orders of magnitude in duration, up to the record observed 
tdur = 1100 S. 

13) The sub-s duration bursts have a values between 20 and 200, 
thus overlapping with the a range of type I burst. Precursor 
events of sub-s durations have been observed before, but they 
are always followed by long type I bursts ^in’t Zand et al.[20l4l (. 
Therefore, the sub-s bursts are of type IT 

14) The relaxation-oscillator behaviour is more accurately pre¬ 
dicted by the burst duration than the fluence (Fig.[^ bottom). 


6 DISCUSSION 

Having laid out all observational features of the type II bursts of the 
RB, we now discuss the physical implications of these. We leave 
the discussion of the instability models for type II bursts to the next 
Section. 


6.1 Accretion versus ejection 

Two observables might suggest the possibility of mass ejection in 
the RB. One comes from the type I burst behaviour, which we dis¬ 
cuss in Sec. [6.4[ The second finding suggesting that mass could be 
ejected during type II bursts is that their peak flux never, for any 
persistent flux or burst fluence, surpasses a certain value which, 
within a factor of 2, is equal to the expected Eddington limit for a 
NS at the distance of the hosting globular cluster, with X = Q.7 and 
an isotropic radiation field (Figs.|^through|^. One might therefore 
wonder whether any surplus radiation power is transformed to the 
kinetic energy of an outflow enforced by radiation pressure. 

There is, however, no spectral evidence for expansion of the 
emitting region. Thus, the picture arises that when the Eddington 
luminosity is reached, the accretion is temporarily halted, inter¬ 
rupting the release of gravitational energy. This in turn lowers the 
luminosity so that accretion can ensue again, rising back to the Ed¬ 
dington limit, and so on. The radiation pressure may momentarily 
cut off the fuel line for the radiation (the accretion flow), but not 
result in isotropic expansion. As a consequence of this, all energy 
is liberated through radiation, and there are no kinetic losses. This 
is in contrast to the type I bursts, where reaching the Eddington lu¬ 
minosity does not affect the energy release from the layer that has 
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been heated by the burning, and the photosphere expands. Instead, 
once a type II hurst flux reaches the Eddington luminosity, a larger 
fluence Eh can only he liberated by means of an increased dura¬ 
tion, which makes for a tight relation between burst duration and 
fluence (Fig.|^ bottom). 


6.2 Two kinds of type II bursts and three kind of modes 

The data point to two kinds of type II bursts: long, often luminous 
ones, and short ones. Long bursts, predominantly mode 0, occur at 
relatively larger mean accretion rates, while shorter ones abound 
at lower rates, during mode 2. They mix in mode 1. Long bursts 
induce humps and adhere more tightly to a relaxation oscillator be¬ 
haviour, while the short ones are not followed by humps and show 
a larger scatter around the relaxation oscillator relation. 

The evidence for this is provided in Fig.[^ At large (Fpers) 
one observes both long and energetic bursts, that become longer 
and stronger at larger persistent fluxes, and short and weak ones, 
displaying an opposite trend (i.e. becoming shorter and weaker at 
larger (T'pers)). The critical fluence separating these two opposite 
feedbacks is Eh ~ 2x10®® erg s“®. This duality is reflected by 
the way in which the bursts affect the persistent emission back. 
While short, weak burst do not have a feedback on Fpers, the 
long and strong bursts that appear in mode 0 and at the end of 
a mode 1 sequence generate humps of enhanced persistent emis¬ 
sion. Fig.fT^shows that the threshold for this behaviour is again at 
Eh ~ 2xL0®® erg s“®. 

Another look at Fig. shows that the threshold in flu¬ 
ence at which the bursts become Eddington limited is Eh ~ 
5x10®® erg s“®. Notice that some mode 2 type II bursts reach 
fpeak = TBdd, but Very quickly decay, and are not followed by 
humps. 

Possibly then, a sustained period of Eddington-limited emis¬ 
sion is necessary for a type II burst to be followed by a period of 
enhanced persistent emission. It would therefore make more sense 
to speak of “shorf’ and “long” bursts, rather than the three modes 
that are mentioned in literature. What Fig.[T^also shows, however, 
is that there is a continuity in behaviour between the two, rather 
than a clear cut distinguishing two different behaviours among type 
II bursts. 


6.3 The relaxation-oscillator behaviour 

Our data confirm the existence of the well-known relaxation os¬ 
cillator between the burst fluence Eh and the time to the next burst 
(Fig. § middle). However, we also find a fighfer correlation 
between the burst duration tdur and At (Fig.[^ bottom). It is quite 
possible that the relaxation-oscillator relation is a by-product of this 
tighter correlation and of the aforementioned one between Eh and 
longer fdur- This would mean that it is the timescales, rather than 
the energetics, that connect the bursts and the persistent emission. 

The burst properties seem to depend weakly on the average 
intra-burst persistent flux (Fig. |^. Also, relatively small changes 
in the persistent flux are visible between modes. The variations in 
the burst properties are therefore much larger than those in the ac¬ 
companying persistent emission. In other words, the variation in 
the global mass accretion rate that must be present between the be¬ 
ginning and the end of the hard state is not equally split between 
the bursts and the persistent emission, but almost entirely goes into 
the former. 
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Figure 14. The recun'ence time tree of type I bursts at different stages in 
the outburst. Blue points are soft-state measurements from B13, from which 
the fit is derived tree <x {dashed line). Squares are averages for 

tree during mode 0, 1 and 2 phases, with the upper and lower limit derived 
using only securely identified type I bursts and including candidates, re¬ 
spectively (with a very small difference during mode 2). They are plotted 
against both the intra-type II burst average flux (Fpers) and (Ftot), which 
includes the burst fluence in the average. The latter does not match the rela¬ 
tion extrapolated from the soft state, but it is more c ompa tible with a steeper 
tree oc FPere relation (dotted-dashed line, see Sec. 6.4 ) 


6.4 Type I bursts as a probe of mass accretion by type II 
bursts 

Our data of type I bursts provide a test of the hypothesis that type 
II bursts are an accretion phenomenon. They provide an indepen¬ 
dent constraint on the accretion rate by comparing type I recurrence 
times with {Epera} and (Etot)- Fig. [l4| shows that the type I 

burst recurrence time in the hard state matches the extrapolated fit 
to the soft state if (Fpera) is taken as a proxy of the mass accretion 
rate, which traces the persistent emission only. Adding the burst flu¬ 
ence to yield (Ftot) gives instead too long recurrence times, even if 
one includes bursts for which the type I identification is tentative. In 
other words, taking into account the mass that type II bursts should 
be accreting on the NS surface does not yield the predicted recur¬ 
rence times. There are three possible explanations for the mismatch 
between Fee and (Etot)- 

Firstly, the type II bursts could be ejection rather than accre¬ 
tion events. The disproportionate coincidence in mode 1 of type I 
bursts with humps during type II burst intermission would then be 
explained by having mass accrete onto the NS only during these 
intermissions and not during type II bursts. Ejection is possible in a 
strong propeller regime, where the disc is truncate d by the magne - 


tosphere at sufficiently large radii (rm > 2®^®rcor, Liietay2014|f 


The mechanism is similar to the one highlighted by [D’An g elo &| 
|Spruit|p010||2012l see below) with matter piling up and accreting 
in burst-like cycles. In this case, accretion is accompanied by the 
launch of a well-collimated, magnetically dominated Poynting jet. 
However, the predicted burst durations and recurrence times fall in 
the ms to s range, and very large B and small mass accretion rates 
(well below 10 per cent of the Eddington rate) are necessary. But 
the strongest argument against this hypothesis is that, in a propeller 
regime, the luminosity should drastically decrease with respect to 
the accreting regime, because the potential energy that would oth¬ 
erwise be liberated during the fall onto the NS surface will not be 


© RAS, MNRAS 00Q,[TpQ| 
















14 Bagnoli et al. 


available jP’ Angelo et al.|2014| >. As the type II bursts reach lumi¬ 
nosities up to the Eddington limit, we are inclined to reject ejection 
as a valid hypothesis. 

Secondly, the type I bursting rate could respond not only to the 
instantaneous mass accretion rate, but also on the rate at which it 
itself varies (i.e. m). Possibly, the sudden and short-lived increase 
in m associated with a type II burst will not produce the same effect 
in the burning layers as if the same mass had been more slowly 
accreted. To the best of our knowledge, however, simulations of 
how the type I burst behaviour changes with accretion rate only 
include models in which m is constant over long periods of time 
(e.g., |Heger et al.|2007| >. It would perhaps be interesting to see what 
these models predict in case of swift, large and short-lived changes 
in the accretion rate, but we cannot presently elaborate further on 
this hypothesis. 

Thirdly, the shorter duration of type I bursts in the hard state 
indicates that, as the mass-accretion rate decreases, the RB ap¬ 
proaches a pure-He burning regime. In this case a steepening of 
the relation in Fig. is expected. Theoretical models predict a 
steeper relation in a pure-He regime ( [Cumming & Bildsten|2000| ), 
compatible with the tree oc measured in IGR J17480-2446 
([Linares et al.|201^ see their fig. 11). We plot this empirical re¬ 


lation in Fig. 14 assuming pure He burning takes place for bursts 


below 10 per cent of the Eddington limit, at the state transition. This 
way we can account for the recurrence times plotted against (Ftot) 
during mode 1 and 2, although type I bursts remain under-abundant 
for mode 0. 

This last hypothesis has the advantage that it also explains the 
observed change in the type I burst duration. As already mentioned, 
type I bursts in the hard state are shorter than those in the soft state, 
both in total duration (~ 20 s, versus ~ 100 s in the soft state) 
and in the duration of their peaks (see Fig.|12|l. This is also ascrib- 
able to changes in the fuel composition: as the mass accretion rate 
decreases, the recurrence time of type I bursts lengthens, leaving 
less H unbumt by the hot CNO process, which in turn reduces the 
burst duration due to the lack of unstable H burning via the slow 
rp-process ( jFujimoto et al.|1981^ . This means that in the hard state 
the RB is approaching a pure-He burst regime. We think this is the 
most likely explanation of the three. 


7 MODELS FOR TYPE II BURSTS 


Since the discovery of type II bursts in the Rapid Burster, an assort¬ 
ment of instability mechanisms have been proposed for their origin. 
Bursts have been variously attributed to an instability in the accre¬ 
tion flow (either in viscosity or temperature), instabilities around 
the innermost stable circular orbit or as a result of the interaction 
between the accretion disc and a stellar magnetic field (e.g., [KuT] 
jpeKj^Ki^njx|£994[ |Spruit & Taam|1993| (. In their review article, 
Lewin et al. ( 1993^ give an overview of the various strengths and 
weaknesses of the different models. We proceed to review these 
models again and test their predictions against our observations of 
the population. 


7.1 Type II bursts: a magnetic phenomenon? 

7.1.1 Theory of disc-magnetosphere interactions 

A weakness of many models proposed to explain type II bursts is 
the apparent uniqueness of the RB, since most proposed instabil¬ 
ities should be widely applicable. In this respect models based on 


magnetospheric instabilities are perhaps the most promising, since 
a strong magnetosphere introduces new characteristic lengthscales 
to the system, and the distinctiveness of the RB behaviour could 
then be attributed to an unusual magnetic field geometry (such as 
alignment with the rotation axis, e.g., |Kuijpers & Kuperus|19^4l l, or 
unusual relationship between the physical properties that determine 
the magnetic regulation of accretion and the properties of the ac¬ 
cretion flow itself (in particular the accretion rate at which the state 
transition occurs). On the other hand, since X-ray pulsations have 
never been detected in the RB, it is not clear that there is a strong 
magnetic field in the system (although by the argument above this 
could be due to a near alignment between the rotation and the mag¬ 
netic axes). The BP does have an estimated magnetic field and 
measured spin period {B ~ (2 — 6) x 10^° G, P = 467 ms, 
jPegenaar et al.|2014[|Finger et al.|1996) . The observation of type 
I X-ray bursts in the RB puts an upper limit of 73 ~ 10^'^ G for 
the field (based on measured fields of know bursters; jPatruno ^ 
|Watts|20T^ , but such a field would still be strong enough to trun¬ 
cate the accretion disc and regulate the gas flow in the inner regions 
of the star at moderate accretion rates (e.g., [Pringle & Rees|1972| (. 
Magnetospheric instability models (studied both analytically and 
via numerical simulations) have also received the most theoretical 
development in the last 20 years (since the [Lewin et al. TT^ re¬ 
view) so we concentrate our discussion on those results and how 
they compare with observation. 

Accretion on to stars with strong magnetic fields becomes reg¬ 
ulated by the magnetic field once it reaches a certain distance from 
the star (the ‘magnetospheric’, or ‘Alfven’ radius, Vm). At the 
disc is truncated by the magnetic field, and the gas is channelled 
along magnetic field lines to accrete on the star jPringle & Rees[ 
[1972[ [Shakura & Sunyaev[[1973f . The magnetospheric radius is 
typically approximated as the point where the ram pressure of the 
infalling gas equals the magnetic pressure of the field (with some 
considerable uncertainty from the fact that both the disc and the star 
are rotating; e.g. [Ghosh & Lamb[1979[[Sptuit & Taam[1993[ l: 


Trrr ~ 3.1 X 10* 
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The location of the inner disc radius is considerably uncertain and 
can vary with the spin rate of the star jSpruit & Taam[[T993[ l. A 
factor ~2 uncertainty for Cm should generally be assumed. 

Since in general the disc and the magnetic field do not rotate at 
the same rate, the coupling between the two will twist the magnetic 
field lines, creating a torque that allows for angular momentum ex¬ 
change between the disc and the star. The sign of the torque will 
depend on the location of Cm relative to the co-rotation radius, the 
radius at which the star’s rotation frequency fl* matches the Kep- 
lerian frequency of the disc: 



( 2 ) 


Accretion can only proceed easily when rin < rc, that is, 
for relatively high accretion rates, slow stellar spin rates, and low 
magnetic field strengths. When rin > Vc the stellar magnetic field 
presents a centrifugal barrier to accretion, and either expels the 
infalling gas in an outflow (the ‘propeller’ scenario; [Illarionov &[ 
[Sunyaev[1975[ l or keeps it confined in the inner disc outside Vc (the 
‘trapped disc’ scenario; Sunyaev & Shakura|1977[[Spruit & Taam[ 


[I^ [D’Angelo & Spruit|2'oT^ 
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7.1.2 Trapped discs 


The trapped-disc picture forms the bases of the most developed 
model to explain type II bursts in the RB. In this scenario, when 
the accretion rate falls below a certain rate, Cm (eq. [TJ moves out¬ 
side Tc, and gas is no longer able to accrete on to the star. As long as 
Vm remains close to Tc, the disc-magnetosphere interaction will not 
provide enough energy to unbind the gas from the disc and launch 
an outflow. Instead, the gas stays confined in the disc, gradually 
piling up as more gas is accreted from larger radii. This situation 
continues until enough gas has piled up to push against the mag¬ 
netic field for Cm to move inside rc. Here the centrifugal barrier 
disappears and the accumulated gas can rapidly accrete until the 
reservoir of gas has been depleted and again moves outside r^. 
Observationally, the instability will manifest as short outbursts (as 
gas accretes freely through the disc and onto the star) contrasted 
with periods of non-accretion (as gas gradually accumulates) with 
a roughly stable luminosity generated by the interaction between 
the disc and the magnetic field. 


This instability was described and studied analytically and nu¬ 
merically in Spruit & Taam| ( |1993) >; |D’ Angelo & Spruit| ( |20ldl l, and 
[D’Angelo & SpruTt[ i 2012^ , who investigated its general properties 
(duty cycle, frequency, outburst shape, etc.) as a function of the 
mean accretion rate and the details of the disc-field interaction. In 
[D’Angelo & Spmit|j20I0| > and [D’ Angelo & Spruit[ ( [20T^ , the disc- 
field interaction is characterized by two (uncertain) lengthscales: 
Ar, the radial extent of the disc that is coupled to the magnetic 
field and Ar 2 , the range of Cm around rc where the disc changes 
from being a trapped disc (outside rc) and accreting freely (inside 
rc). A key feature of this instability is that it produces bursts on a 
wide range of timescales, which can generally be much longer than 
the dynamical timescale of the inner disc, and plausibly match the 
burst timescales seen in the RB. 


In the trapped disc instability, the main driver of changes in 
burst shape and duration is a changing mean accretion rate into the 
inner disc. Roughly speaking, the recurrence time increases and 
the burst duration and fluence decrease as the mean accretion rate 
decreases. This is broadly seen in the RB: the longest bursts in 
mode 0 occur early in the outburst, and the time between mode 
2 bursts increases gradually as the overall luminosity declines. If 
this instability is responsible for the behaviour seen in the RB, the 
stochastic behaviour seen in the variation between individual bursts 
is mainly attributable to a variable accretion rate in the inner disc 
regions. This may be consistent with the X-ray spectral transition 
from a disc-dominated (‘soft’ or ‘high’) state to a power-law domi¬ 
nated (‘intermediate’ to ‘low’) state. In neutron star and black hole 
binaries, the spectral transition to a power-law spectrum is typi¬ 
cally accompanied by large luminosity fluctuations on a range of 
timescales, which are commonly interpreted as fluctuations in lo¬ 
cal mass accretion rate propagated through the accretion flow (e.g., 
[Lyubarslai][ 1991) . Since the instability requires an accumulated 
mass reservoir, these accretion rate fluctuations will manifest as a 
variation in the duration and frequency of bursts. 

However, the simulations of [Spruit & Taam[ ( [1993^ and 
[D’Angelo & Sprult| j2010[ l do not show the characteristic ‘relax¬ 
ation oscillator’ behaviour seen in type II bursts in the RB, where 
the wait time between bursts scales with the fluence of the previous 
burst. This is because neither set of simulations properly accounted 
for viscosity in the large-scale disc, which sets the timescale to re¬ 
fill the inner disc regions. Instead the simulations supplied gas into 
the inner regions of the disc at a steady rate, which leads to a steady 
burst period for a given accretion rate. In a realistic disc, if the reser¬ 


voir of gas is emptied before the outer disc can resupply it with gas, 
the build-up time will be longer than has been observed in simula¬ 
tions. This introduces a new timescale into the problem: the viscous 
refilling time. If the accretion rate into the inner disc is unsteady, 
the timescale between bursts will be variable, but will be at least 
as long as it takes for gas to flow into the inner disc. This may 
account for the recurrence time-fluence relationship seen in Fig.[^ 
the recurrence time has a minimum (the viscous inflow rate) but can 
be much longer (if the supply rate suddenly decreases). [D’Angelo[ 
[& Spruit] ( [2012[ l indeed found that a steadily decreasing accretion 
rate (in a large disc) can induce chaotic accretion bursts from the 
trapped disc instability (fig. 10 of that paper). This hypothesis re¬ 
mains to be investigated further. 

[D’Angelo & Sprult|p012| > identified two distinct instability 
regions for the trapped disk instability: ‘RT (long duration, low 
amplitude bursts), and ‘RIF (shorter, smaller amplitude bursts). In 
‘RF bursts, Ar 2 /r < 0.02 (i.e. the disc makes abrupt transitions 
between the accreting and non-accreting states), and the bursts of¬ 
ten have complex profiles and with strong contrasts between burst¬ 
ing and inter-burst luminosity. In contrast, the ‘RIF bursts occur 
at larger values of Ar 2 , bursts are more sinusoidal in shape with 
higher luminosities, shorter periods and less contrast between burst 
and inter-burst luminosities. The complex burst profiles and char¬ 
acteristic outburst timescales of the mode 0 bursts suggest the RI 
instability is more likely applicable for the Rapid Burster. 

At high mean accretion rates, the ‘RF accretion bursts are 
longer, with burst profiles that frequently show quasi-periodic os¬ 
cillations (see Fig. [T^ at the tail of the outburst, similar to mode 
0 bursts of the RB. A flat-topped burst is also sometimes seen, al¬ 
though typically with an initial spike of accretion, which is not seen 
in RB mode 0 bursts. As the mean accretion rate drops, the burst 
duration and fluence decreases, and the bursts become more widely 
separated in time, similar to evolution in the mode 2 bursts of the 
RB (see Fig.[T](. 


7.1.3 Application to the RB 

In comparing the predictions of the trapped disc model to the RB 
bursts, we consider only the short burst (mode 2 and short mode 1 
bursts). This is because the long mode 0 bursts have significantly 
different properties: they are often Eddington-limited bursts with a 
distinctive burst profile and persistent emission that shows a hump 
and (occasionally) type I bursts, suggesting significant accretion 
onto the star continues between bursts. We speculate that for mode 
0 bursts, the trapped disc instability could appear in conjunction 
with a second physical process (perhaps related to the dramatic 
state change that marks the onset of type II bursts), but we do not 
consider this further. 

Applying the trapped disc burst model to the RB is challeng¬ 
ing, as its spin rate and magnetic field are unknown. However, if 
we assume (as this model would suggest) that the type II bursts 
begin when Vm moves just outside rc, then the ‘critical accretion 
rate’ (that is, for the onset of bursts) is about 0.1 MEdd. Equat¬ 
ing Cm and Tc (equations!^ andat this accretion rate leads to a 
degeneracy between B and P*. We can set some bounds on these 
parameters from other evidence: the appearance of type I bursts 
suggests B < 10^° G, while P* > 0.0017 s to make the RB com¬ 
patible with the fastest-rotating known neutron stars. This leads to 
a range of P ~ 6 x 10^ — 10^° G for spin periods P ~ 2 — 80 ms 
(with the fastest spin periods corresponding to the smallest mag¬ 
netic fields). The smaller rc, the nearer the disc will extend to the 
star, and the shorter the recurrence time will be between bursts. 
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Figure 15. Bursts predicted by the trapped disc instability for a neu¬ 
tron star with B = 10^® G and a spin period P = 0.02 s, for dif¬ 


ferent mean accretion rates (see Sec. 7.1 . From top to bottom, M = 


[0.03,0.005,0.002] MBdd, with the burst fluence decreasing {F = 
[7.8,1.9,1.8] X 10®* erg) and recumence time increasing as the mean ac¬ 
cretion rate decreases. We plotted the properties of these simulated bursts 
against those of our sample in Fig. |7] The viscous response of the disc to 
the changing inner disc may create a delay in refilling the inner disc for 
larger bursts, reproducing the behaviour seen in the RB. Note that between 
bursts accretion in the inner disc is suppressed, and the disc luminosity is 
generated by the interaction between the disc and the magnetic field. 


because the burst timescale is generally set by the viscous diffu¬ 
sion timescale, fvisc oc r^^/v in the disc’s i nner regions. For an 
a-disc viscosity, v ~ a{H/r)^{GMr)^^^ i Shakura & Sunyaev 


|1973| where H is the scale-height of the disc), so that for a ~ 0.1 
(H/r) ~ 0.1, fvisc ~ 10-1000 s at for the range of possible 
spin periods. In |D’ Angelo & Spruit| ( |201^ , the strongest bursts oc¬ 
cur on timescales 0.1-10 fvisc, which would correspond to 1-100 s 
if the RB is a rapid rotator and 100-10"^ s if the spin frequency is 
lower. There is thus considerable degeneracy in the model, but for 
a reasonable range of spin periods, magnetic fields and viscosities, 
it can reproduce the observable features of the RB. 

In Fig. [T^ we show bursts at three different accretion rates, 
choosing representative simulated parameters that reproduce the 
RB range of luminosity within the constraints on spin and mag¬ 
netic field set above. We adopt B = 10^° G, P = 0.02 s 
and a viscosity v = 10“®(GMr)“®, which gives strong bursts 
when the mean accretion rate decreases below ~ 10 per cent 
Lndd- From top to bottom, the accretion rate decreases: M = 
[0.03, 0.005,0.002] MEdd- The flux properties of the model agree 
quite well with the observed properties of the RB (see Fig.|^. The 
fluence in the bursts decrease from 7.8 x 10®® erg in the brightest 
burst to 1.8 X 10®® erg in the weaker bursts. In comparison, the RB 
shows bursts over a somewhat larger range of fluence, with 90 per 
cent of bursts between 0.13 -1.8 x 10®® erg. The peak flux in the 
simulated bursts is nearly independent of accretion rate, which is 
also roughly consistent with the RB, where the peak luminosity of 
90 per cent of the bursts spans a range of ~ 2. 

As discussed in Sec. |4.1[ the persistent emission between 
bursts does not show obvious correlations with other burst prop¬ 
erties, and varies by a factor ~ 2 for the majority of bursts. This is 


naturally accounted for by the trapped disc model, where between 
bursts accretion in the inner disc is suppressed, and the disc lu¬ 
minosity is generated by the interaction between the disc and the 
magnetic field, which adds angular momentum and energy to the 
disc. Since the location of the inner disc does not vary much be¬ 
tween bursts, the luminosity also stays roughly constant, regardless 
of the fluence of the subsequent burst. As mentioned above, this 
does not account for the ‘humps’ and occasional type I bursts seen 
between long mode 0 bursts, since the distinctive nature of mode 
0 bursts suggests additional physics beyond the theoretical picture 
outlined above. 

The observed range of burst recurrence times and durations 
are less well accounted for in the model. As described above, all 
simulations of the instability assumed a constant accretion rate in 
to a small region of the inner disc, which does not properly account 
for the time it could take the disc to refill after an outburst. As a 
result, the recurrence time scales linearly with accretion rate - the 
lower the mean accretion rate the longer it takes to refill the reser¬ 
voir and the weaker the subsequent burst. This is the opposite of 
what is apparently observed in the RB. However, if the accretion 
rate in the inner disc is strongly variable (as is postulated to explain 
rapid luminosity variability in the hard state; |Lyubarskii|1997| l, the 
relationship between the mean accretion rate and burst fluence is 
less well-defined. Simulations also show bursts with a smaller vari¬ 
ation in burst durations (about a factor of 4, with durations between 
~ (2 — 8) s) than the RB, in which the majority of bursts (excluding 
mode 0, which as we argue above likely have an additional physi¬ 
cal process governing their behaviour) last between 2-24 s, with a 
peak around 3 s. 

In summary, the trapped disc instability can reasonably ac¬ 
count for multiple properties of the RB short type II bursts for a 
variable accretion rate in the inner disc regions. Simulations can 
match observed burst fluences, maximum luminosities, and persis¬ 
tent flux properties for reasonable (albeit poorly constrained) spin 
and magnetic field. The range of simulated burst durations and flu¬ 
ences are somewhat smaller (by a factor ~ 4 and ~ 2 respectively) 
than is observed in the RB. Most significantly, the relaxation os¬ 
cillator behaviour is not reproduced, although as we argue above, 
this is at least in part because the response of the large-scale vis¬ 
cous disc to the burst instability has not been properly modelled. 
The finite viscosity in the disc will delay the refilling time once 
gas is depleted, lengthening the time between bright bursts. The 
RB uniqueness may be explained by a special combination of mag¬ 
netic field strength and stellar spin period (which together sets the 
accretion rate where the instability can manifest) and an alignment 
between the magnetic field and the spin axis. 

7.2 Other magnetospheric instabilities 

Several other mechanisms have been proposed to produce instabili¬ 
ties in the accretion flow that could explain the RB behaviour, but it 
is not clear that the timescales on which they are expected to oper¬ 
ate will match the RB. Both mechanisms outlined below rely on the 
ideal MHD properties of the accretion flow, i.e. that field lines are 
frozen in to the plasma and may not efficiently re-couple to the disc 
once the connection is severed (which is expected to happen rapidly 
due to differential rotation between the field and the disc). This can 
lead to either compression of the magnetic field by the disc mov¬ 
ing inwards or unstable fingers of accretion onto the neutron star 
surface via interchange instabilities. 

[Kuijpers & Kuperus| jl994^ proposed that the perfectly con¬ 
ducting disc can be completely shielded from the magnetic field of 
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the star, so that the disc will spiral inwards and pinch the dipolar 
magnetic field along the equator into a distorted ‘butterfly’ shape. 
This compression lasts until the matter touches the surface of the 
star, whereupon the gas penetrates through the field and accretes on 
to the star. This releases the pressure on the magnetic field, which 
rapidly re-expands and restarts the cycle. Some version of this sce¬ 
nario is seen in all magnetic accretion simulations (e.gJMiller^ 
Stone|1997||Goodson et al.|1997||Hayashi et al.|1996||Romanova 

et al.|2009> , which use ideal MHD and hence shield the disc from 

the magnetic field. However, due to the finite numerical diffusivity 
present in all simulations, the timescale for the instability is typi¬ 
cally a few dynamical timescales, which is much shorter than the 
RB bursts (~ 50 — 1000 dynamical times, depending on the spin 
of the star). Moreover, it is not at all clear theoretically that gas will 
have difficulty coupling to the magnetic field, as is required by this 
picture. Simulations of MRI turbulent discs with stratified density 
profiles (e.g. [Davis et al.|2010^ show the disc field extends verti¬ 
cally out of the disc, which could form a low-density magnetized 
corona above the disc surface (e.g. [Jiang et al.|2014| ). This would 
provide an obvious site for reconnection, which could recouple the 
disc to the star and prevent strong compression of the magnetic 
field. 

In this scenario, the burst should have a gradual rise (as the 
accretion disc spirals inwards toward the star) with a sudden spike 
when the gas finally hits the surface. This also seems inconsistent 
at least with the short bursts, which show low persistent emission 
between bursts, although it could perhaps explain the inter-burst 
humps seen between mode 0 bursts. In order to match the observed 
bursting timescales, the original [Kuijpers & Kuperu7| ( [1994[ l paper 
also posited a strong 10^^ G field for the RB, which is inconsis¬ 
tent with the observation of type I X-ray bursts at relatively high 
luminosities. 

An alternative scenario will develop when the field is not 
strong enough to efficiently disrupt the disc, allowing accretion to 
proceed via interchange (Rayleigh-Taylor) instabilities on to the 
star. In this picture the RB magnetic field could be significantly 
weaker than assumed above (or the spin rate lower), so that the disc 
extends down to the star until the accretion rate drops to ~ 10 per 
cent I/Bdd- At lower luminosities the field will truncate the disc 
but not enough to enforce accretion on to the poles, and fingers of 
accretion will extend down to the stellar surface. This was studied 
in depth in the in the 3-D simulations of Kulkami & Romanova 
( [Kulkarni & Romanova|2008l l. However, as in the above scenario, 
the timescale for the instability is typically the orbital timescale of 
the inner disc edge, which will be smaller than the spin period of 
the star and thus likely considerably shorter than the observed burst 
timescales unless the RB is a very slow rotator. 


7.3 Thermal/viscous instability models 

A number of models have also been proposed that do not rely on 
a role for the magnetic field. [Lewin et al.| jl993[ > have already re¬ 
viewed them and explained their shortcomings; we will briefly dis¬ 
cuss them in the light of our sample properties. 

A class of models invokes thermal ([S hakura & Sunyaev|1976^ 
or viscous instabilities ( [Lightman & Eardley[1974[ ) in the accretion 
disc to explain the luminosity fluctuations. Both instabilities arise 
in the inner portion of an a disc, where the radiation pressure domi¬ 
nates and electron scattering is the main source of opacity. Thermal 
instabilities are due to an inability to maintain a thermal balance, 
while viscous instabilities arise because the viscous stress becomes 


inversely proportional to the surface density, breaking up the inner 
disc in concentric rings as a result. 

These analytical results however only established local cri¬ 
teria for instabilities to grow. They were first tested in a global, 
time-dependent simulation by [Taam & Lin[ ( [T984[ ), who confirmed 
the existence of unstable solutions. The produced bursts liberate 
about the right amount of energy (10®®“^® erg), and the observed 
durations (0.1-1 s) and recurrence times (~ 10 s) at least overlap 
with the (however much broader) ranges observed in type II bursts. 
Also, the instability develops in a range of mass accretion rates 
((0.1 — 0.45) X 10^® g s“^, between 10 and 45 per cent I/Edd), 
above which the flow becomes stable again. 

There are two main problems with this model. Firstly, differ¬ 
ent time scales and energetics can only be reproduced with varia¬ 
tions of the global mass accretion rate. Although this problem is 
shared by the class of magnetospheric models discussed earlier, we 
should stress that at least in that case a much broader variety of 
burst properties is observed, closer to matching the ranges of flu- 
ences, durations and recurrence times that we have characterized in 
this study. Secondly, models relying on disc instabilities leave the 
uniqueness of the RB unaddressed, as they should apply to all ac¬ 
creting compact objects, even black hole binaries, and have indeed 
been put forward to explain the variability of sources like Cyg X-1. 


7.4 GR models 


An entirely different approach was taken by [Walk^ ( [1992^ , who 
proposed that the RB be an exceptional source in hosting a NS with 
a radius smaller than the innermost stable circular orbit (ISCO), 
surrounded by a massive accretion disc. The latter gives the very 
low viscosity the model requires. Also, the magnetic field must be 
weak enough not to disrupt the accretion flow. In this case, rather 
than having the magnetosphere act as a gate, radiation torques in¬ 
side the ISCO lead to a sudden spike in the accretion rate, that re¬ 
peats after the inner disc has been replenished on a viscous time 
scale. 

The predicted burst luminosities match those observed in type 
II bursts, and the bursts last between roughly 1 and 50 s. A very 
large disc would further be compatible with the large orbital period 
of the BP (Porb = 11.8 d), while no orbit information is avail¬ 
able for the RB. The model suffers from the same problem that 
afflicts all other models, namely the lack of a relaxation oscilla¬ 
tor behaviour. But its main shortcoming is that it cannot account 
for type II bursts in the BP, which was only discovered three years 
after this model was proposed. |Walker[ ( [1992) claim that the pres¬ 
ence of a massive disc means that even relatively large magnetic 
fields (up to 10^^ G) would not disrupt the Keplerian disc before it 
crosses the ISCO. Although estimates of the magnetic field strength 
in the BP are as low as (2 — 6) x 10^° G i Degenaar et al.|2014 1 , 
there cannot be any doubt that magnetospheric accretion is hap¬ 
pening in that source, because of the presence of pulsations at the 
spin frequency, both in the persistent emission and in the bursts. If 
GR effects were truly responsible for type II bursts in the RB, then 
another mechanism altogether would be at play in the BP. 


8 CONCLUSIONS 

Our population study has highlighted a number of new, impor¬ 
tant features of the type II bursts. The bursts are Eddington lim¬ 
ited, which leaves little doubt as to their origin lying in accretion 
episodes. We have constrained the duration of these instabilities to 
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cover the 10“^ — 10® s range, and we have shown that a mini¬ 
mum recurrence time of order of 10 s seems necessary to develop 
a burst. We encourage dropping the old mode 0, 1 and 2 classifica¬ 
tion, adopting instead a simpler one between short and long bursts, 
which reflects the dichotomy these show in their interaction with 
the persistent emission. 

We judge models based on thermal and viscous instabilities 
too general to explain the scarcity of sources displaying type II 
bursts, while models based on GR effects imply that a different phe¬ 
nomenon altogether is responsible for the BP bursts. Models based 
on the gating role of the magnetosphere can reproduce a number 
of properties in our sample, while at the same time requiring a 
fine tuning between the stellar spin, the magnetic field strength, the 
mass accretion rate and the alignment between the magnetic field 
and the spin axis, and thus account better for the uniqueness of the 
RB. 
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APPENDIX A: BURST SEARCH ALGORITHM 

The thousands of type II bursts from the RB in the RXTE archive necessitate 
an automatic burst seai'ch algorithm for the identification process. We now 
proceed to explain the algorithm that we developed from trial and error. 
During the development, results were verified by visual inspection of the 
light curves. 

The automatic procedure employs a light curve generated at 1 s resolu¬ 
tion from all PC A ’Standard-I’ data with the RB in the FOV (2.4 Msec), 
for all photon energies and for PCU2 which is the only PCU that is al¬ 
ways on. Furthermore, the pointing and the elevation of the source above 
the Earth’s horizon was tracked. From the pointing, the off-axis angle be¬ 
tween the optical axis to the RB and 4U 1728-34 were calculated. Pointing 
and elevation information is available at 16-s resolution. For unsampled 
times, the pointing information was taken from the prior data point. 

First, for each light-curve bin an estimate was made of the persistent 
emission underlying a potential burst at that bin, by determining the average 
of all flux values within 100 s prior to that bin, provided that the pointing is 
constant within 0.05 deg for at least the previous 5 s. If there was no reliable 
background from those data, the average flux was taken from the posterior 
20 to 100 s interval. The background thus determined was subtracted from 
the flux of each bin. Second, a ’significance’ was determined by taking the 
ratio to the Poisson error of that bin. If the significance was larger than 
1, the bin was further investigated for the presence of a burst. Third, the 
background-subtracted flux of subsequent bins after that bin were added to 
the signal and the accrued significance determined. This process was re¬ 
peated for as long as the significance increased, allowing for it to decrease 
only once. When finished, an attempt was made to accrue signal going back¬ 
ward in time from the initial bin. Fourth, the candidate burst was verified by 
checking that no other bursts previously determined (i.e., starting at earlier 
times) were in the new burst interval, that the total significance was at least 
10, that the elevation was at least 2.5 degrees, that the pointing was con¬ 
stant during the burst and that there was at least one ’empty’ bin between 
two consecutive bursts. A list of candidate bursts was thus composed. Since 
this allows for a more careful background determination from the original 
flux an'ay by excluding the burst intervals and, fifth and last, the procedure 
was repeated. When it was impossible to determine the background, the 
burst was ignored. 

All type I bursts from the RB show the typical fast-rise-exponential- 
decay (FRED) shape. Soft-state type I bursts have rise times (7.1 s) and 
durations (~100 s; B13), which are relatively long for type I bursts, and 
typical of the burning of H-rich fuel ^Lewin et al.|1993) . As we have shown 
in Sec |4.3| hard-state type I bursts are typically shorter (~40 s). Nonethe¬ 
less, type II bursts of comparable durations cannot mimic the FRED shape, 
because multiple peaks appear in their decays jTan et al.|1991^ . This allows 
identification of a RB type I burst even when its profile in the light curve 
overlaps with those of several fast-recurring type II bursts (as in Fig |13| 
bottom). 

As for type I bursts from 4U 1728-34, they show instead the much 
faster time scales (rise time <1 s, duration ~10 s) that are typical of H- 
poor nuclear fuel jGalloway et al.||2Q10^ , and based only on this could 
be confused with short RB type II bursts. However, the bursts from 
4U 1728-34 show a characteristic bimodal distribution of peak fluxes 
(9.2x 10~® erg cm~^ s~^ for the bursts showing photospheric radius ex- 
pansion (PRE) and 4.5x10“^ erg cm“^ s“^ for normal bursts, [Gallo way | 
|et al.|2Q()^ , which are much larger than observed in RB type II bursts. Still, 
when the PCA is pointed at the RB the peaks fluxes of the weakest 4U 1728- 
34 bursts can be similar to those of the brightest RB type II bursts. Again, 
the multiple peaks in the decays reveal the nature of the latter. 

Thus, a list of 8754 bursts was derived. Eighteen of these concern 
overlap between two kinds of bursts. The breakdown of three kinds of bursts 
is: type I bursts from the RB (123, 9 overlapping), type-I bursts from 4U 
1728-34 (167, 9 overlapping) and type II bursts from the RB (8458, 18 
overlapping). Twenty-two bursts could not be definitely identified because 
of insufficient statistics. They are either type I from the Rapid Burster or 
type II. Otherwise, type-I bursts from 4U 1728-34 are easily recognizable 
because they are approximately two times brighter than the brightest type 
II bursts from the RB and last the shortest. The latter characteristic is nec¬ 


essary for identification because, when the PCA points at the RB and 4U 
1728-34 is an off-axis angle of 0.5 degrees, the collimator shadow cuts the 
measured peak flux in roughly two and is comparable to the brightest bursts 
of the RB. For the other two pointing configurations, identifying bursts from 
4U 1728-34 is unambiguous. The difference between type I and II bursts 
from the RB is also manageable. Type I bursts are at least 20 s long, have a 
typical long tail jBagnoli et al.|2Qf3]|2014) and show no ringing, contrary 
to most type II bursts. 

The start and end times of the bursts were determined by searching, 
starting from the bin with the highest flux, downward and upward, respec¬ 
tively, for the last times when the flux drops below 10% of the peak flux. 
Notably, we identified a record duration for a mode-0 burst. Observation ID 
92026-01-06-03 apparently starts off during a type II bursts which subsides 
only after 1100 s in a fashion typical for a mode-0 burst. 


APPENDIX B: SPECTRAL ANALYSIS 

In principle one would like to perform a full time-resolved spectroscopic 
analysis of each type II burst to determine parameters such as peak flux 
(in erg cm~^s~^), burst fluence (in erg cm“^), fluence of the persistent 
emission between bursts and a. However, given previously published spec¬ 
tral analyses of RB data d ata (e.g. , [Marshal l et al.|1979| [Rutledge et al.| 
1 1995 1 [Marshall et al.|2Q0l]|Mahasena et al.|2Q03) , there is merit in asking 

whether that is necessaiy. Spectral variability during and between type II 
bursts is presumably limited and the average energy per photon is probably 
nearly constant. In order to verify this and determine a conversion factor 
between PCA-detected photons and energy, we analyzed a representative 
subset of the data. The subset was drawn from the ’offset’ observations that 
are free from contamination by 4U 1728-34. Eleven observations, identi¬ 
fied by unique ObsIDs, were chosen representing typical as well as extreme 
behavior of the RB. They ai'e listed in Table [Bl| 

The time resolution of the spectral analysis varied across the different 
observations and within observations is chosen to follow the variability in 
flux. In total 374 spectra were extracted from these 11 observations, involv¬ 
ing 8 mode-0 bursts (including the longest one detected), 25 mode-1 bursts 
and 23 mode-2 bursts. 

The spectral analysis involves the following generalities. Spectra were 
extracted from event-mode data with 64 spectral channels. All active PCUs 
were employed. The bandpass was limited to 2.5-20 keV which is where 
the PCA is well calibrated jjahoda et al.|200^[Shaposhnikov et al.|2012) 
and contains most of the signal. A systematic eiTor of 0.5 per cent per chan¬ 
nel was included (following [Shaposhnikov et al.|2012) . All spectra were 
corrected for particle background and cosmic diffuse background as pre¬ 
dicted through the PCABACKEST tool (version 3.8). The PCA response ma¬ 
trix was calculated with version 1 1.7.1 of tool PCARSP. This includes a cor¬ 
rection for the collimator response which is always non-standard because 
the RB is always considerably off-axis in the selected ObsIDs. The collima¬ 
tor throughput for the RB was usually about 0.4. Spectral bins were com¬ 
bined to obtain at least 20 counts per bin to ensure applicability of the 
statistic. All spectral modelling was carried out with XSPEC version 12.8.2. 
A constant cold interstellar absorption component was as sumed to apply, 
equivalent to Wh = 1-6 x 10^^ c m~^ (Frogel et al.|l995f . The absorption 
model of jMomson & McCammQn|jl983) was followed. 

A significant contribution may be expected from the Galactic Ridge. 
To determine this, off-axis observation ID 92026-01-04-00 was chosen 
which represents the observation with the lowest flux measured for the RB 
field. It was taken on June 9, 2006, 11 days before the onset of the out¬ 
burst presented in Fig. 1. The background-subtracted data was modelled 
with the combination of a Raymond-Smith component jRaymond & Smith 

B and a power law, following the prescription of [Valinia & Marshall 
as also determined from PCA ridge data. The fit was satisfactory 
with xt ~ 1-79 (i' = 20). This constant spectral component was assumed 
to be present in all studied RB spectra. 

We tested various models on the data: single component models 
(power law, Planck function, thermal bremsstrahlung, Comptonization) and 
models with combinations of two components. It was found that mod¬ 
els with three free parameters were sufficient to describe all spectra. We 
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Table Bl. A representative subset of RB obsei'vations for spectral studies 


ObsID 

Date of 

Observation 

Type of behavior 

Bolometric flux range 
(10“^ erg cm“^s“^) 

204I8-01-01-01R 

17-jun-1997 

bright soft state 

1.2 

20093-01-07-01R 

26-jun-1997 

high soft state 

1.3 

20418-01-07-00 

lO-jul-1997 

fainter soft state 

0.26 

30424-01-01-02R 

25-aug-1998 

mode 0, short bursts 

0.1-1.7 

30424-01-02-02R 

29-aug-1998 

mode 1, sequence of 5 bursts 

0.1-2.5 

30424-01-03-02R 

l-sep-1998 

mode 1, sequence of 20 bursts 

0.1-2.8 

30424-01-04-02R 

4-sep-1998 

mode 2, sequence of 22 bursts 

0.1-0.8 

30424-01-06-01R 

lO-sep-1998 

mode-2, single burst in tail outburst 

0.1-0.8 

92026-01-06-00 

23-jun-2006 

mode 0, sequence of 5 short bursts (brightest fluxes from RB) 

0.2-3.3 

92026-01-06-03 

26-jun-2006 

mode 0, record-long burst (>1300 s) 

0.15-0.8 

92026-01-06-04 

27-jun-2006 

temporary return to soft state 

0.45-0.55 
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Figure Bl. Time-resolved spectroscopy for ObsID 92026-01-06-00. The 
response is corrected for the off-axis position of the RB. The fitted model is 
the Comptonization model according to|Titai'chukH1994). 


chose the Comptonization model according to |Titarchuk| j 1 99^ ; |Titarchuk| 
|& Lyubarskij|(T995) as baseline since this is the simplest model fitting all 
data. Furthermore, we find that this model performs best to the highest qual¬ 
ity spectrum (20093-01-07-01R whose spectrum contains 10® counts), with 
_ Q 92 for 1 / = 20. As an example, for the same data x? = 11-3 for a 
model consisting of a disk black body and a power law, and 6.88 for a cut¬ 
off power law. We note that it was necessary to complement the continuum 
model with a narrow line at about 6.5 keV, in the spectral region where the 
Fe-K resonance line resides (6.4-6.9 keV, depending on the ionization state 
of Fe). Possibly it is related to instrumentally scattered emission from the 
bright source 4U 1728-34 just outside the field of view, or to an inadequate 
modelling of the Galactic Ridge at this location. 

Figure [BTj shows an example of a time history for observation 92026- 
01-06-00 which shows a sequence of 5 bright and relatively short mode-0 
bursts. What is clearly visible here, and that applies in general to all our 
spectral results, is that plasma temperature and optical depth show no sig¬ 
nificant variability, while photon seed temperature does. The plasma tem¬ 
perature varies in general between 2.5 and 3.5 keV while at the same time 
the plasma optical depth varies between 7 and 4, both typical 1-sigma errors 
of 1 in value. A more detailed analysis of the spectral data is subject for a 
future paper. We here confine the analysis to infening an average translation 
factor between photon fluxes and fluences and energy fluxes and fluences. 

For the most accurately determined 112 spectral shapes, we find that a 
photon flux of 1 c s“^PCU“^ for the whole bandpass translates on average 
to a 2.5-20 keV energy flux of (1.02 di 0.02) x 10“^^ erg cm“^s“^. All 
112 conversion values are contained within 7% from this average. Thus, 
the maximum systematic error by assuming a constant conversion factor is 
1%. Converting from 2.5-20 keV energy flux to unabsorbed bolometric flux 
involves a factor of 1.24± 0.07. All 112 bolometric correction factors range 
between -7% and -1-2% from the mean. Therefore, the maximum systematic 
error in the bolometric correction factor is 7%. 

In Fig. |B2| we show spectra of the five most diverse states of the 
RB. The conversion factor from c s~^PCU~^ to unabsorbed bolometric 
erg cm~^s~^ varies between 1.23 x 10“^^ and 1.56 x 10“^^. 
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Energy (keV) 

Figure B2. Unfolded 2.5-30 keV v — spectra for five widely different 
states of RB. From bottom to top: non-burst emission between slow mode-2 
bursts (black curve, ObsID 30424-01-06-01R), non-burst emission in long 
intermission between mode-1 bursts (light blue, 30424-01-02-02R), peak 
of longest mode-0 burst ever measured (green; 92026-01-06-03), soft state 
prior to mode-0 (red; 20418-01-01-01R), peak of bright short mode-0 burst 
(dark blue; 30424-01-01-02R). The lowest spectrum shows the iron line 
clearly whose origin is undetermined. 
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